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Abstract 

The profile of Navarro, Frenk, and White (the NFW profile), which was derived from the TV-body 

m 

^v^j ■ simulations of cold dark-matter halos, is a strong candidate for a galaxy or cluster profile. In order to check 

t— i ; 

ON . the usability of the NFW profile as a first approximation of a galaxy model, we studied the characteristic 
ON 

overdensity and scale radius of galaxies by reproducing the image positions and flux ratios of 2-image 
gravitational lens systems, under the following simple assumptions: the galaxies are spherically symmetric, 
and stars and external shear do not contribute to the gravitational lens. The scale radii of the lcnsing 
galaxies are smaller, and the characteristic overdensities are larger than the predicted value in the TV-body 
simulations. These results indicate that our assumptions are overly simplified. It may be impossible to 
simply adopt the NFW profile, which does not include stars, to probe the cosmological parameters or the 
light propagation in an inhomogeneous universe and so on. If we adopt a softened isothermal profile to the 
lensing galaxies, the scale radii and the central matter densities agree with models which are derived from 
other observational results for early-type galaxies and which are independent of gravitational lensing. The 
isothermal sphere as a first approximation of a galaxy model has no serious problem. 
Key words: dark matter — galaxies: structure — gravitational lensing 
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1. Introduction 

Gravitational-lens systems provide a powerful tool to probe the potentials of galaxies, the cosmological parameters, 
light propagation in an inhomogeneous universe and so on. In an analysis of gravitational-lens systems or other 
observations, we need to assume a profile of the mass distribution of galaxies. The softened isothermal profile and the 
modified King profile have been conventionally used to determine the Hubble constant, H , or the density profiles of 
lensing galaxies (Kent, Falco 1988; Grogin, Narayan 1996; Keeton, Kochanek 1997). These profiles are approximately 
derived from the collisionless Boltzmann equation, and explain the H I rotation curves well (Binncy, Trcmaine 1987). 
If these profiles, however, do not represent the actual features of galaxies, the cosmological parameters obtained need 
to be altered. Thus, it is important that we correctly select the mass profiles of galaxies. 

There is no definite evidence that an isothermal profile is not appropriate. We should, however, compare the 
isothermal profile with other plausible profiles. One of the most highly possible candidates for galaxy profiles, 
instead of the isothermal profile, is the profile described by Navarro, Frenk and White (1995). Using iV-body 
simulations of cold dark-matter halos, Navarro, Frenk and White (1996, hereafter NFW96; 1997, hereafter NFW97) 
revealed that a density cusp forms at the center, and that the density profiles have the same shape, independent of 
the halo mass and of the initial density fluctuation spectrum. According to Navarro et al., the central density profile 
fits the analytic form, 



(hereafter we refer to this as the NFW profile), where r s , p cr = 3Hq/8ttG, and 5 C are the scale radius, the critical 
density and the characteristic overdensity, respectively. The NFW profile is found in 7V-body/hydrodynamic simu- 
lations (Navarro et al. 1995; Eke et al. 1997). They studied objects with masses ranging from those of dwarf galaxy 
halos to those of rich galaxy clusters. Carlberg et al. (1997) found from measurements of galaxy velocity dispersions 
in clusters that the NFW profile can reproduce the measured velocity dispersion. 

The results for the scale radius and the characteristic overdensity are summarized in table 1 in NFW96. Navarro 
et al. predict that the scale radius is approximately several hundred kpc at the cluster mass scale. On the other 
hand, X-ray observations of galaxy clusters do not agree with these results (Makino et al. 1998; Makino, Asano 
1999). The size of r s for galaxy clusters predicted in NFW96 is smaller than the observed values. The relation 



p(r) = 



(r/r s )(l + r/r s ) 



2' 



(1) 
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between the scale radius and the characteristic overdensity of galaxies also needs to be probed. 

In general, a large number of stars or gas should be incorporated when we obtain accurate parameters of lens 
galaxies, although the NFW profile is obtained from cold dark-matter simulations. Complex models with many 
parameters can greatly reproduce the observations of gravitational lenses. However, it is not absolutely clear that 
the rough parameter values of the dark halo can be obtained from the NFW profile, which does not include stars or 
gas. If we can use the NFW dark halo as a first approximation of a galaxy model, the NFW profile would be very 
convenient to apply it to various cosmological analyses. On the other hand, the isothermal profile has been solely 
adopted as a first approximation. The light propagation in an inhomogeneous universe have studied by applying an 
isothermal profile to galaxies or galaxy clusters (e.g. Jaroszyhski 1992; Premadi et al. 1998, Tomita et al. 1999). In 
this kind of analysis, we need simple models of the galaxy profile, while a complex model with many parameters and 
perturbations is needed to obtain the Hubble constant from gravitational lenses. The parameter values of the NFW 
profile has been predicted by numerical simulations. If we can neglect the contribution of stars, light propagation 
in an universe with galaxies of the NFW profile can be promptly studied. An analysis with the NFW profile may 
bring us different results from those at present for light propagation. For the above reasons, the usability of the 
NFW profile, which does not include stars as a first approximation of the galaxy profile, should be checked, as is 
the case with the isothermal sphere. The usability of the NFW galaxy halo as a cosmological tool is not necessarily 
well-known. In an analysis of galaxy clusters with the X-ray surface brightness (Makino et al. 1998), the gravity 
of the gas is not very important in order to obtain the density profile. Of course, this situation is not the same for 
galaxies, since their mass-to-light ratio is much smaller than those of galaxy clusters. We expect that the mass-to- 
light ratio in the central regions of galaxies is about 10 (e.g. Lauer 1985). Sofue (1999) claimed that the bulge of 
LMC has a large fraction of dark matter (the mass-to-light ratio is about 20-50). For spiral galaxies, a high-density 
concentration of dark mass was reported (Sofue et al. 1997, 1998). Thus, it is not necessarily obvious that the stellar 
component predominates in the centeral region of galaxies sufficiently enough to intervene the simple analysis of the 
dark halo profile. We should confirm whether the NFW halo, without stars, can be a rough approximation of galaxy 
profile or not. 

In this paper we consider the characteristic overdensity and the scale radius of galaxies for the NFW profile from 
2-image gravitational lens systems and a simple lens model without stellar component. The softened isothermal 
profile is also adopted for a comparison. This profile represents the other conventionally used profiles. The reason 



4 [Vol. , 

for studying this subject is to confirm whether the NFW profile, which does not include stars or gas, can be 

a cosmological tool to a first approximation. We discuss whether the NFW profile, which is obtained from our 

method, can express the actual order of galaxy parameters or not. Our goal in the future is to check the consistency 

between observations and the parameters predicted from the JV-body simulations. This study is the first step toward 

this goal. Only the predicted parameter values in NFW96 are constraints for the parameters of the NFW profile at 

present. Bartclmann (1996) and Evans and Wilkinson (1998) adopted a scale radius which satisfies an estimation by 

numerical simulations of NFW96 for galaxy clusters. Thus, our strategy in this paper is that we accept the results 

in NFW96 to be as Bartclmann et al. assumed, and to compare the parameters obtained from the gravitational 

lenses with those in NFW96. The disagreement between NFW96 and the results of gravitational lenses indicates 

that the galaxy models which do not include stars are not adequate, or that the results in NFW96 are not reliable. 

If the results in NFW96 are wrong due to softening and small numbers of particles, it also means that we cannot 

use the NFW profile for galaxies as a first approximation at present. Anyhow, we can clarify the usability of the 

NFW profile from this study. We apply simple lens models, with no contribution of stars, no external shear, and 

spherically symmetric lensing, as a first approximation. Our aim is not to obtain accurate parameters of galaxies by 

a precision analysis. The usability of the NFW dark halo as a first approximation is our theme in this article. In 

section 2 we explain the method of analysis. We only consider spherical lenses so as to avoid complex assumptions. 

2-image lens systems are good targets for this method (Schneider et al. 1992). In section 3, our method is applied to 

the observation of lens systems. The results for the NFW and the softened isothermal profiles are summarized. The 

results for the softened isothermal profile are compared with the galaxy models of Jaroszyhski (1992). In section 4, 

we summarize our work, and discuss our results. 

Throughout this paper, we assume H = 50 kms" 1 Mpc" 1 , f^o = 1, and A = 0. 

2. Method of analysis 

We consider 2-image lens systems in order to estimate the relation between the galaxy parameters. Since 2-image 
lens systems are approximately expressed by a spherically symmetric lens (Schneider et al. 1992), wc can study the 
structure of lensing galaxies by the following simple process without using the \ 2 minimization technique. We do 
not need to consider any complicated assumptions such as ellipticity, external shear, or spiral structure. 
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The spherically symmetric lenses are completely described by their surface mass density, S(r), and the critical 

surface mass density, 

c 2 D 

S cr =- =2-, (2) 

cr 4nG D\D\ s y ' 

where D\ is the angular diameter distance to the lens, D s is that of the source, and D\ s is the distance from the lens 

to the source. By symmetry, the lens position, the source position, y, and the image positions, r, are on one straight 

line. By placing the lens at the origin, the lens equation reduces to a one-dimensional form, and is expressed by 

y = r-—-, (3) 

where m(r) is defined by 

m (r) = 2 f Ep-r'dr', (4) 

JO ^cr 

(Schneider et al. 1992). For the NFW profile, one finds 

m(r) = r 2 K g(x), (5) 
where 



g(x) = 21n | + { 



4 1 x — l , 

r% — 7 tan \— T ^ >:L 
Vx 2 - 1 \ x + l 

1 (a = 1) ( 6 ) 



4 , i /l — a; . 

:tanh J (0 < x < 1), 



I Vl-x 2 V 1 + x 

and x = r/r s (Bartelmann 1996; Maoz et al. 1997). The dimensionless parameter kq is written by 

k q = 2<5 c p cr r s /S cr . (7) 

Here, we neglect the contributions of stars and gas as a first approximation. We also consider the softened isothermal 
profile, 

P(r) = = , ^ V2 ■ (8) 
1 + (r/r s y 

Here, we used the same notation for the parameters as those for the NFW profile, in order to simplify our formulation. 
We should note that there is no direct relation between the parameters for the NFW and the softened isothermal 
profiles. In this notation, m(r) for the softened isothermal profile is the same as equation (ph with 



<J 



(x) =ir(y / l + x 2 -1). (9) 
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In spherical lens systems, two images are positioned on opposite sides of the lensing galaxy. The observed image 

positions > and < should satisfy 
m(r A ) 



'"A 



m{\r B \ 

Irul 



''A 

from equation (||). 

The amplification factor is written by 



(10) 



1 



m(r) 



1 



d m(r) 



dr r 

The observed flux density ratio, /, should satisfy 
j = mOa) 



(11) 



(12) 



From equations (|l0| ) and (|l2|), we can uniquely determine the parameters r s and Ko- The characteristic overdensity, 
5 C , can be obtained from these parameters. If the two parameters, <5 C and r s , are determined, the mass of the halo 
can also be obtained. However, the total mass of the NFW profile is divergent at r = oo. In NFW96, r2oo, which 
determines the mass of the halo, was defined as M200 = 200p cr (47r/3)r| 00 . Namely, the mean density of the halo 
within r2oo is 200p cr . The characteristic overdensity for the NFW profile is analytically expressed as 
200 c 3 



(13) 



3 ln(l + c)-c/(l + c) 
where c = r 2 ooAs- 

For the softened isothermal profile, we conform the method of mass estimation to the method in NFW96. We 
obtain 

^ = -z (14) 

3 c — tan c 

in a similar way. Both S c 's in equations jl^ ) and (|l^) are monotone increasing functions of c. We will discuss the 
parameter c, instead of 5 C . Since M200 is nothing but a criterion in each profile, we do not directly compare the 
masses for the two profiles. This circumstance also holds for the other parameters also. 

In NFW96, the characteristic overdensity of the halos with masses M200 was estimated to range from those of 
dwarf galaxy halos to those of rich galaxy clusters. The results are summarized in table 1 in NFW96. We compare 
our results with those of NFW96. In addition, we calculate the time delay, the image position of the undetected 
third image C and its flux, using the determined parameters r s and d c . 
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3. Results 
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3.1. Samples 

As a result of our simple assumptions for lensing galaxies, the number of usable lens systems is limited in our 
estimation. Few 2-image lens systems have been discovered. In addition, the lensing galaxies of the lens systems 
have not all been found. We exclude the double quasar 0957+561 from our samples because a cluster of galaxies 
apparently contributes to the wide separation between the images of this system (Young et al. 1981). The observed 
parameters of our samples are listed in table 1. 

We adopt the data of flux density ratio in the reddest or longest wave band in order to minimize the influence 
of dust extinction. The flux-density ratios are less well determined than the image positions, because of the source 
variability, microlensing, extinction, and so on. To account for this, we broadened the flux error bars to 20% in a 
similar method to that of Keeton and Kochanek (1997). The lens and image positions are not exactly on a straight 
line in the observational data of the 2-image lens systems. Generally, the observational error of a lens position is 
larger than one of the image positions. Therefore, in order to adopt our method, we force the lens position to be 
placed on the nearest point to the observed lens position on a straight line joining the two images. We obtained the 
image positions in a one-dimensional form, and list them in table 1. The distances from the observed lens positions 
to the moved lens positions are 0".002, 0".035, 0".121 and 0".08 for B1600+434, B1030+074, SBS 1520+530 and 
HE 1104—1805, respectively. Some of them are larger than the observational errors of the lens positions. This fact 
indicates that we need external shears in our analysis of these systems. 

When an external shear is weak, it slightly moves the image positions, and breaks the alignment between the 
images and lens. We can assume that the error of the image positions in the spherical approximation, due to the 
external shear, is on the same order of the lens-position shift in the weak external shear limit. In this case, the errors 
of the obtained parameters due to the error of the image positions are smaller than those due to the flux error, as we 
discuss in the next subsection. Thus, when estimating the order of the parameters of a galaxy profile, the effect of 
external shear is negligible. We neglect the errors of image positions, which is equivalent to neglecting the external 
shear as a first-order approximation. In this method, the fitness of the statistic cannot be evaluated. Alternatively, 
taking into account the flux error, we estimate the uncertainties in the results. 

The lensing galaxy of B1600+434 is an edge-on spiral (Jaunsen, Hjorth 1997), and that of B1030+074 has a 
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substructure (Xanthopoulos et al. 1998). We, however, neglect these structures in our calculation. 



3.2. Mass Profile of Galaxies 

The results of our calculation are summarized in table 2. The errors of the results in table 2 are derived from the 
flux error bar. We neglect the errors of the image positions due to external shear. The shift of the lens position in 
SBS 1520+530 is the largest one at 0".121. The errors of the image and lens positions in this method of analysis are 
on the same order as the shift of the lens position. We force the lens position in SBS 1520+530 to shift 0".l from 
the position given in table 1 toward images A or B, and calculate the parameters of the lens galaxy for the NFW 
profile. Next, we move the both positions of images A and B, in table 1, 0.05 arcsec (a total of 0".l) toward or away 
from the lens position, and calculate the parameters. In these cases, the maximum differences from the results in 
table 2a in the scale radius and the characteristic overdensity are 0.33 kpc and 0.20 x 10 7 , respectively. These are 
smaller than ones due to the flux error in table 2a. The effects of the errors due to the image positions of the other 
lens systems are much smaller, and negligible. This situation is the same for the isothermal profile. Thus, we can 
neglect the errors of the image positions in this method of analysis. 

HE 1104—1805 has no solution for the NFW profile in equations (10) and (12). This means that our methods 
cannot explain this lens system. The system is an extraordinary system. The mass of the lensing object is large, 
because of the high redshift of the lensed images and the wide angular separation between the images. The redshift 
of the lensing object is very high at z — 1.6616 (Courbin et al. 1998). A reason for the absence of a solution might be 
that the system is influenced by a cluster of galaxies or external shear. This possibility remains to be proved because 
no obvious overdensity of galaxies is detected in the immediate region surrounding the lens (Courbin et al. 1998). 
This system has a solution for the softened isothermal profile. However, the predicted flux of the undetected third 
image C is brighter than that of image B. Therefore, our simple analysis would be undesirable for HE 1104—1805. 
HE 1104—1805 needs some perturbations. 

On the other hand, the softened isothermal profile can not explain the image positions and flux density ratio of 
the lens system B1030+074, because of the extremely high flux ratio. The flux-density ratio is one of the highest 
for any known lensed system (Xanthopoulos et al. 1998). However, the NFW profile can explain this lens system. 
For B1030+074, the NFW profile is more favorable than the isothermal profile, unless one takes into account the 
substructure of the lens galaxy. 
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The two profiles of the mass distribution have their merits and demerits for each observation of gravitational lens 
system. We cannot determine which profile is superior from only the gravitational lens. As a result, we conclude 
that the reproducibilities of image positions and flux ratios are comparable for the two profiles. 

In figure 1, we plot the dimensionless parameter c against the scale radius r s for the NFW profile in tabic 2a (filled 
square). The results of numerical experiments of M 2 oo < 10 14 M Q in NFW96 are also plotted (open square). Using 
X 2 minimization, they are empirically approximated by 



for M 200 < 1O 14 M . This equation is a nothing but a gauge to make easier to compare our results with those of 
the simulations. Apparently, the values of parameter c of the three lensing galaxies are larger than the results of 
NFW96, and the scale radii, r s , of the lensing galaxies are smaller than those predicted by NFW96. The large value 
of c means a large characteristic overdensity. The largest value of c in figure 1 is derived from B1030+074, the lens 
system of the highest flux density ratio. We also plot c against the mass of the halo, M 20 o, in figure 2. The relation 
of these two parameters in NFW96 is approximately expressed as 



for M 2 oo < 10 1 Mq. Figure 2 is equivalent to figure 1, because the NFW profile is a 2-parameter model of the 
matter distribution, figure 2, however, clearly indicates that the values of c arc much larger than those of NFW96 
for the same value of M 2 oo- We have compared the parameters derived from the simulations at z = with those of 
the high redshift lensing galaxies. Since the central concentration of galaxies becomes larger, in general, as galaxies 
evolve, the high-redshift lensing galaxies should have smaller values of c than the predicted value at z = 0, as is in 
NFW96 or NFW97. Therefore, the result, that c of the lensing galaxies is too large, is unchanged. 

For the softened isothermal profile, a diagram for c versus r s in table 2b is shown in figure 3. The lowest value of 
c in figure 3 is the estimate for HE 1104—1805. The core radius is very large. We also plot for the galaxy models of 
the softened isothermal profile described in Jaroszyfiski (1992). These models have often been used to study light 
propagation in an inhomogeneous universe (Jaroszyhski 1992; Premadi et al. 1998). The characteristic overdensity 
in the models of Jaroszyhski is given by 




(15) 




(16) 
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where v is the velocity dispersion. The velocity dispersions are given by 




1/4 



km s 



(ellipticals) , 



1/4 



V 



= < 



km s 



l 



(SOs), 



(18) 



0.381 



km s 



(spirals), 



where L and are the luminosity and characteristic luminosity, respectively. The Faber-Jackson relation (Faber, 
Jackson 1976) and the Tully-Fisher relation (Tully, Fisher 1977) were used for early-type and spiral galaxies, respec- 
tively. The dependence of the scale radius on the luminosity has been studied by Lauer (1985), Kormendy (1987), 
and Fukugita and Turner (1991). Their studies are based on surface photometry profiles, H I rotation curves, or 
the virial theorem which are independent of gravitational lensing. In conformity to these studies, the models of 
Jaroszyhski adopt 



In figure 3, all data points of the gravitational lenses, including HE 1104—1805, are obviously close to the lines 
for the models of Jaroszyhski for early-type galaxies. Although the lensing galaxy of B1600+434 is a spiral, the 
models for early-type galaxies are favorable. The value of c given by the model for spirals is smaller than those of 
the lensing galaxies. r s of HE 1104—1805 is too large. This system needs some perturbations as we have mentioned 
above. However, the orders of c in the other two systems are available. On balance, there is no serious problem 
using the isothermal sphere as a first approximation of a cosmological tool. 

We cannot directly compare the consistencies of the NFW and the softened isothermal profile. However, the 
difference between the results of NFW96 and ours is conspicuous, while the models of Jaroszyhski well account for 
the parameters of the lensing galaxies. As we shown above, the softened isothermal sphere as a first approximation 
of a cosmological tool has no serious problem in an analysis of light propagation in a universe. We cannot use the 
NFW profile, without stars, and need more complicated assumptions in such analysis. 

The predicted time delays of B1600+434 and SBS 1520+530 for the NFW profile in table 2a are about 10% larger 
than those for the isothermal profile in table 2b. This different time delay for the two profiles will sufficiently affect 
the estimate of the Hubble constant. 




2 I — I kpc, (spirals). 



200 — pc, (ellipticals and SOs), 



(19) 



Eliminating L/L* by equations ( |is| ) and (f[9|), the c-r s relation can be plotted. 
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4. Conclusions and discussion 
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We applied the NFW and the softened isothermal profile to lensing galaxies. We assumed that there is no contri- 
bution of stars and external shear, and that the lens galaxies are spherically symmetric, as a first approximation. By 
reproducing the image positions and their flux ratios, we estimated the scale radii and the characteristic overdensi- 
ties. The reproducibilities of image positions and flux ratios are comparable for the two profiles. If we assume that 
the mass distribution of galaxies is described by the NFW profile, the scale radii are smaller, and the characteristic 
overdensities are larger than those in NFW96. If we adopt the softened isothermal profile to the lensing galaxies, the 
scale radii and the central matter densities are not so different from those in the models of Jaroszyhski for early-type 
galaxies, which are derived from the various observational results. The central densities of the Jaroszyfiski model 
for spiral galaxies are too small to account for the observations of the lensing galaxies, including the spiral lens in 
B1600+434. The predicted time delays for the NFW profile are substantially larger than those for the isothermal 
profile. 

If we uncritically accept these results, they are more advantageous to the softened isothermal profile than to the 
NFW profile. However, there are two possible and natural explanations of our results for the NFW profile. One is 
that our lens model may be overly simplified. We assumed that the contributions of stars and gas are negligible. 
A large amount of stars may influence the image deflection. We cannot simply adopt the NFW profile, which does 
not include stars, to probe the cosmological parameters or the light propagation in an inhomogeneous universe. We 
need more complicated models with stars in their analyses. We also assumed that the external shear is weak. If the 
external shear is strong, our method cannot be applied. 

Another explanation is that the results in NFW96 for galaxies are wrong due to softening and small numbers of 
particles. In this case, much higher resolution is required in A-body simulations in order to describe the parameters 
of galaxies. For the present, we cannot decide which explanation is more accurate. Therefore, we need much 
more sophisticated lens models which include the contribution of stars or external shear, which we plan to develop 
in the future. Anyhow, the NFW profile for galaxies, which does not include stars, cannot be applied as a first 
approximation of galaxy model at this stage. The isothermal sphere as a first approximation of a galaxy model has 
no serious problem. 
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We discuss other possible explanations for our results below. In NFW96, the smallest gravitational softening is 

1.5 kpc. This size could be too large to estimate the scale radius of galaxies. However, the small gravitational 

softening may influence the central profile of galaxies in TV-body simulations. The power of the density cusp at the 

center is controversial because of effects of gravitational softening and poor resolution at small scales. Fukushige 

and Makino (1997) showed from simulations that as particle numbers increase 1 order of magnitude, typical central 

density profiles become shallower than p oc r~ 2 , but steeper than p oc r _1 . Moore et al. (1998) also indicate that 

the central density profile of cold dark-matter halos becomes steeper as the resolution in the cosmological TV-body 

simulations becomes higher. The density profiles of Moore et al., or other steep profiles, need to be adopted to the 

lcnsing galaxies. These are our future objectives. 

As mentioned in subsection 3.2, the effects of galactic evolution cannot remedy the contradiction of the character- 
istic overdensity, because the central concentration of galaxies generally becomes larger as the galaxies evolve. The 
high-redshift lensing galaxies should have smaller values of c than the predicted value. 

If we assume that the initial density fluctuation spectra, P{k) oc k n , is not restricted by <r%, the rms mass 
fluctuation in spheres of radius 8/i _1 Mpc, we can enlarge the value of c in TV-body simulations. While the standard 
model, er 8 = 0.63, was assumed in NFW96, a universe with various power spectra in NFW97, n = 0, —0.5, —1, 
and —1.5, produces larger values of c. As index n becomes larger, the value of c becomes larger. From figure 9 in 
NFW97, universes with indices n = 0, or —0.5 seem favorable for our results on gravitational lensing. 

If the cosmological constant, A , has a finite value, we can marginally decrease the value of c for the lensing galaxies. 
For example, c's values are 77.7, 299, and 68.1 for B1600+434, B1030+074, and SBS 1520+530, respectively, in the 
flat A universe; H = 75 kms -1 Mpc -1 , £! = 0.25, and A = 0.75. Here, we expressed A in units of 3Hq. However, 
judging from figure 9 in NFW97, these results do not agree with the TV-body simulations. 

Of course, the statistical reliability of our estimation is not high enough, because of the small number of samples. 
As mentioned, the number of 2-image lens systems is not very large. More sophisticated lens models can be adopted 
to more complex lens systems and we will have to increase the number of samples. 

We are grateful to Nobuyoshi Makino and Kenji Tomita for their helpful advice. This work was supported by a 
Research Fellowship of the Japan Society for the Promotion of Science. 
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Figure Captions 

Fig. 1. Diagram for the dimensionless parameter c against the scale radius r B for the NFW profile. The filled squares are 
derived from our results in table 2a. The error bars derived from the 20% error of the flux ratio are also plotted. The 
open squares are the results of the Af-body simulations for M200 < 10 14 Mq in table 1 in NFW96. The solid line is the 
best-fit relation between c and r B for the results of NFW96. 

Fig. 2. Diagram for the dimensionless parameter c against the mass M 2 oo for the NFW profile. The filled squares, the open 
squares, and error bars are the same as in figure 1. The solid line is the best-fit relation between c and M200 for the 
results of NFW96. 

Fig. 3. Diagram for the dimensionless parameter c against the scale radius r s for the softened isothermal profile. The filled 
squares were derived from our results in table 2b. The error bars are the same as in figure 1. The solid line, dotted 
line, and short dashed line indicate the models of Jaroszynski for elliptical, SO, and spiral galaxies, respectively. 
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Table 1. List of 2-image gravitational lens systems.* 



Name 


21 


Z s 


TA 




/ 


References 








arcsec (kpc) 


arcsec (kpc) 






B1600+434 


0.4144 


1.589 


1.14 (7.5) 


-0.25 (1.6) 


1.25 


1,2 


B1030+074 


0.599 


1.535 


1.37 (10.4) 


-0.16 (1.3) 


14 


1,3 


SBS 1520+530 


0.765 


1.855 


1.19 (9.7) 


-0.38 (3.1) 


1.89 


4 


HE 1104-1805 


1.6616 


2.316 


1.14 (9.6) 


-2.06 (17.4) 


3.47 


5 



* z\ and z B are redshifts of the lens and the source, respectively, references.- (1) Fassnacht and Cohen (1998); (2) Koopmans, 
de Bruyn, and Jackson (1998); (3) Xanthopoulos et al. (1998); (4) Crampton, Schechter, and Beuzit (1998); (5) Courbin, 
Lidman, and Magain (1998). 



Table 2(a). Results for the NFW profile.* 



Name 


r s 


5c 


C 


M200 


rc 


/CB 


Time Delay 




kpc 


10 7 




10 12 M© 


arcsec 




days 


B1600+434 


9 oc+0.34 
z,.OO_ 28 


1 4Q+034 
^•^-0.30 


92.6± 7 ;2 


i - U8 -0.09 


-o.i6i° ;°oi 


40+0 07 


59.2 + _2'.o 


B1030+074 


58+ ' 12 
u -°°-o.io 


69.9^;° 


372±l\ 


u.oyo_ 055 


-0.04±g;gi 


n7+ - 05 


207.2+2. 7 


SBS 1520+530 


o QK+1.37 
°-° d -0.90 


i 01 +0.95 
iiOi -0 .59 


88.311?;? 


9 01 +0.65 


-o.iii : o3 


U.10_Q Q 5 


1^1. 1 _ 12 6 


HE 1104-1805 








no solution 









* rc is the image position of the undetected third image, /cb is the flux density ratio of image C to B. The time delays of 
the image B lagging behind A are also listed. 
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Table 2(b). Results for the softened isothermal profile. 



Name 


r s 


5c 


C 


M 20 o 




fcB 


Time Delay 




pc 


10 7 




10 13 M© 


arcsec 




days 


B1600+434 


51911° 


9.77ll;S 




1 46 +0 ' 04 


-0 17+ - 01 

U- 1 ' -0.02 


u.ou_ og 


49.3l°; 7 7 


B1030+074 








no solution 








SBS f520+530 


71 c+254 
' 1J -232 


o 9C+8.99 




9 Qfi+ ' 23 


u - lz -0.04 


n 15+° 10 


H7.5«; 4 


HE 1104-1805 


1 Q2no +440 ° 

lazuu -3900 


D 17Q+0-057 
a -0 .038 


163^ 


183±|i 


77+ - 04 

u - ' ' -0.06 


9 40+0.68 


9/19 9+65.6 
88.8 



i i i r~n 




j i i i i i i i i i i i i i i i i i i i i i i i i 

1 10 1CT 

Scale radius [kpc] 



